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We present two examples of velocity distributions for light dark matter particles that reconcile 
the annual modulation signal observed by DAMA with all other negative results from dark matter 
searches. They are: (1) a conventional Maxwellian distribution for particle masses of 6 to 9 GeV; 
(2) a dark matter stream coming from the general direction of Galactic rotation (not the Sagittarius 
stream). Our idea is based on having a signal in Na, instead of I, in DAMA, and can be tested in 
the immediate future by CDMS-II (using Si) and CRESST-II (using O). 

PACS numbers: 14.60.St, 98.80.Cq 



The nature of dark matter is one of the fundamen- 
tal problems of physics and cosmology. Popular candi- 
dates for dark matter are weakly interacting massive par- 
ticles (WIMPs) . Direct searches for dark matter WIMPs 
aim at detecting the scattering of WIMPs off of nuclei in 
a low-background detector. These experiments measure 
the energy of the recoiling nucleus, and are sensitive to a 
signal above a detector-dependent energy threshold. 

One such experiment, the DAMA collaboration pj, has 
found an annual modulation in its data compatible with 
the signal expected from dark matter particles bound to 
our galactic halo 0. Other such experiments, such as 
CDMS 00, EDELWEISS @, and CRESST 0,13, have 
not found any signal from WIMPs. It has been difficult 
to reconcile a WIMP signal in DAMA with the other 
negative results. 

Here we address this issue by posing the following 
question: is there a possible dark matter signal above 
threshold for DAMA and under threshold for CDMS and 
EDELWEISS, so that the positive and negative detection 
results would be compatible? The answer is yes. 

We assume that dark matter consists of WIMPs whose 
mass and cross section we choose in a purely phenomeno- 
logical way, without any attempt of providing an elemen- 
tary particle model to support them. 

The minimum dark matter particle velocity required 
to produce a certain nuclear recoil energy E is 
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where /i = mMj (m + M) is the reduced WIMP-nucleus 
mass, m is the WIMP mass and M is the nucleus mass. 
The nuclear energy threshold E t h observable with a par- 
ticular nucleus corresponds to a minimum observable 
WIMP velocity, the velocity threshold fth- Threshold 
energies for several direct detection experiments are col- 
lected in Table 1. 

To make our argument as simple as possible, consider 
for a moment the case m <C M. Then /1 ~ m is inde- 
pendent of the nucleus mass M, and Vth is proportional 



to 1/ MEth ■ Using the nuclear masses of Na and Ge, 
M Na = 21.41 GeV and M Gc = 67.64 GeV, and the en- 
ergy thresholds in Table 1, the product ME t h is smaller 
for Na in DAMA than for Ge in CDMS-I. For m < M, 
the vth of CDMS-I is 2.44 times that of DAMA, and the 
vth of CDMS-II and EDELWEISS (as well as those of 
other experiments using heavier nuclei) are larger. Using 
the full Eq. J3>, it is easy to see that the velocity thresh- 
old of Na in DAMA is smaller than that of Ge in CDMS-I 
for m < 22.3 GeV. 

This means that for light enough WIMPs it could 
be possible to have a virialized or non-virialized (a 
stream) halo component which could have a velocity 
above threshold for Na in DAMA, and below threshold 
for Ge in CDMS and EDELWEISS. That is, we could 
have a dark matter signal visible for DAMA but not ob- 
servable by CDMS and EDELWEISS. 

CDMS has a small component of Si too. Si is lighter 
than Ge, although heavier than Na, Ms; = 26.16 GeV. 
Given the mentioned nuclear energy recoil thresholds, the 
velocity threshold of Si in CDMS-I is smaller than that 
of Na in DAMA for all WIMP mass values. However, 
considering the CDMS efficiency close to 5 keV energies 
is about 10%, the exposure of the Si detector of CDMS 
near threshold is about 0.3 kg-day, which may be too 
small to have detected the signal which DAMA might 
have seen though its Na. In any event, CDMS has not 
yet used its Si component to set limits on dark matter, 
but only to help in background rejection. 

Light nuclei arc used by CRESST, in particular O, 
M = 14.90 GeV. CRESST-I used sapphire (A1 2 3 ), 
which besides O contains Al, similar in mass to Si. 
CRESST-I has set limits on dark matter with a very low 
nuclear recoil threshold of 0.6 keV, but with a small ex- 
posure of only 1.5 kg-day. The velocity threshold for O 
in CRESST-I is so low that CRESST-I is sensitive to 
the bulk of the halo dark matter particles we are propos- 
ing. CRESST-II uses calcium tungstate (CaW04), which 
also contains the light O nucleus, but background dis- 
crimination sets a relatively high threshold of ~ 10 keV. 
CRESST-II has run a prototype without neutron shield 
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FIG. 1: Stream heliocentric speed u s tr (here assumed to ar- 
rive from ecliptic latitude /3 s tr = 30° and ecliptic longitude 
Astr = 340°) vs. WIMP mass m. The hatched region is com- 
patible with DAMA, CDMS, EDELWEISS, and CRESST. 
Which experiment limits the compatible region is indicated 
along its edges. Also marked are the speed above which 
the stream is extragalactic (dashed horizontal line), and the 
threshold speeds w t h of several experiments and target nuclei. 



and set the limits quoted in Table 1 pj. The completed 
CRESST-II will test our idea. 

Fig. 1 shows the relevant velocity thresholds i>thj as 
functions of the WIMP mass m. Our idea is that a con- 
tribution to the density of WIMPs with velocities smaller 
than the CDMS threshold but larger that the DAMA 
threshold could explain the data. We are here thinking 
in terms of a non-virialized stream of dark matter |9j 
added to an isothermal halo model. The stream would 
be the dominant signal in DAMA (because there would 
be few halo WIMPs above the DAMA threshold) while it 
would change the signal due to the bulk of the halo (e.g. 
in CRESST) very little. The vertical axis of Fig. 1 gives 
the average dark matter stream velocity relative to the 
Sun. Assuming a stream as described below, the region 
where our idea works is the hatched region in the figure. 
To our surprise, the allowed region reaches zero stream 
velocity in a small interval of WIMP masses near 6 to 
9 GeV. In this interval, a stream is not necessary and 
just the contribution of an isothermal halo would do. 

If there is no halo contribution above threshold in 
DAMA, the stream arrival direction is limited by the re- 
quirement that the DAMA modulation peaks May 21 ± 
22 days. May 21 is 61 days after the Spring equinox 
(March 21), and thus the Sun is at ecliptic longitude 
61/365.25 x 360 = 60°. Since the radius vector to the 
Sun and the velocity of the Earth are perpendicular (for 
a circular orbit), the Earth is moving toward a point 
of ecliptic longitude 60° - 90° = 330°. Thus, in case 



the DAMA modulation is entirely due to a dark matter 
stream, the arrival direction of the stream on Earth must 
have ecliptic longitude A = 330° ±22°. The amplitude of 
the modulation depends on the projection of the stream 
velocity onto the ecliptic, and is proportional to cos/3 st r, 
where /3 s t r is the ecliptic latitude of the stream arrival 
direction. 

The non-virialized dark matter stream we are propos- 
ing could be bound to our galaxy or not. There are 
non-virialized dark matter streams bound to our own 
galaxy, such as the tidal streams of the Sagittarius dwarf 
galaxy 10j. The Sagittarius leading tidal stream passes 
through the solar neighborhood, with a heliocentric speed 
of w 350 km/s, but its arrival direction (A, 0) w (187°, 8°) 
does not have the correct ecliptic longitude to give the 
observed phase of the DAMA modulation without a sub- 
stantial contribution from the usual halo component. So 
the stream we are implying should be a different, yet 
undiscovered, stream. There may be dark matter bound 
not to our galaxy but to our Local Group of galaxies > 
and also dark matter bound to our supercluster, possi- 
bly passing through us 12]. Its galactocentric incoming 
velocity v- m is increased by gravitational focusing while 
falling into our galaxy to become a galactocentric ve- 
locity w local = yj v fn + v csc near tne Sun. Here v csc is 
the local escape velocity from the Galaxy. The den- 
sity of an incoming stream is also increased by focus- 
ing, at least linearly with the ratio v\ oca i/vi n , but pos- 
sibly by much larger factors, which are however com- 
plicated to evaluate. Interestingly, the Sun moves rela- 
tive to the Local Group at 306 ± 18 km/s in direction 
(A,j3) = (4° ±10°, 57° ±4°) H3, the ecliptic longitude of 
which is not far from the desired value. 

Our procedure is the following. For the WIMP veloc- 
ity distribution, we take /(v, t) = /h(v, t) + / s tr(v, t), the 
sum of a halo distribution /h(v, t) and an optional contri- 
bution from a dark matter stream / str (v,t). For /h(v, t) 
we take the conventional truncated Maxwellian used in 
the comparison of direct detection experiments, 
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Here = erf(z) — 27r~ 1 / 2 ze~ z , with z = 
normalization factor, v Q is the velocity of the Sun rel- 
ative to the Galaxy (galactocentric velocity), and v®(i) 
is the velocity of the Earth relative to the Sun. For v Q , 
we take the conventional value Vq = 232 km/s, which is 
within the measured value of 233±3 km/s in the direction 
(A,j3) = (341° ±1°, 60.5° ±0.5°) 13]. For v®, we assume 
a magnitude of 29.8 km/s along a circular orbit on a plane 
inclined by 60° with respect to v©. Furthermore, we take 
the conventional velocity dispersion Vh — 220 km/s and 
the escape speed from the galaxy v csc = 650 km/s. We 
make no claim that this is a realistic velocity distribution, 
but it offers us a definite benchmark for comparison. In 
addition, we assume the conventional value for the local 
dark matter density p = 0.3 GeV/cm 3 . With the halo 
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Constraint 


xtex. 


CDMS-I 


Si: 2.83 
Ge: 28.3 


5 


# < 10keV- 7 ft 
_B < 20keV: 22.8 
E > 20keV: 38 


5-55keV: <2.3 events (f) 


[31 


CDMS-II 


Si: 5.26 
Ge: 52.6 


10 


E < 20keV: 22.8 
E > 20keV: 38 


10-100keV: <2.3 events (f) 


r ^1 

[4] 


EDELWEISS 


Ge: 8.2 


20 


100 


20-100keV: <2.3 events (f) 


[51 


CRESST-I 


A1 2 3 : 1.51 


0.6 


100 


(t) 


[7] 


CRESST-TT 


GaWO,,- 10 448 


10 


100 


Ca+O, 15-40keV: <6 events 
W, 12-40keV: <2.3 events (f) 


[81 


DAMA/NaI-96 


Nal: 4123.2 


I: 22 (o) 
Na: 6.7 (o) 


100 


1- 2keVee: <1.4/kg-day-keVee (★) 

2- 3keVee: <0.4/kg-day-keVee (*) 


[61 


DAMA/NaI-03 


Nal: 107731 


I: 22 (0) 
Na: 6.7 (o) 


100 


2-4keVee: 0.0233±0.0047/kg-day-keVee (•) 





TABLE I: Experimental constraints used in this study. Notes to the table: (f) upper limit assuming no detected event; {%) to 
reproduce the published curve in Q, we impose appropriate upper limits all along the recoil spectrum in their Fig. 1; (o) from 
an electron equivalent threshold of 2 keVee, using the quenching factors Q = E ee /E equal to 0.09 for I and 0.3 for Na [l|; (*) 
approximations that reproduce the published a p vs. m limit across our mass range; (•) amplitude of annual modulation. 



model we assume, the maximum possible heliocentric ve- 
locity of a halo particle is v csc + Vq — 882 km/s. 
For the stream, we assume 

/ S tr(v,<) = — e -|v-v str +v ffi (t)|V^r ; (3) 

where £ s tr is the local dark matter density in the stream 
(in units of 0.3 GeV/cm 3 ), v str is the stream velocity rel- 
ative to the Sun (heliocentric velocity), and v s t r is the 
velocity dispersion in the stream (which we fix at 20 
km/s; our results depend very little on the value of Tj str ). 
To produce Figs. 1 and 2, we take £ str = 0.03 and a 
stream arrival direction (A str , /3 str ) = (340°, 30°), so that 
the stream is at 30° of the ecliptic and at 30° of the Sun's 
galactocentric velocity. Such a stream is extragalactic if 
its heliocentric velocity exceeds 840 km/s (marked by the 
dashed horizontal line in Fig. 1). Due to gravitational fo- 
cusing, to obtain wi OC ai > 660 km/s, for example, i.e. a 
stream heliocentric radial velocity v str > 892 km/s, re- 
quires v m > 114 km/s, a reasonable value. 

For each choice of WIMP velocity distribution, we 
must make sure that we not only produce the correct 
amplitude for the DAMA modulation in the viable region 
but also that all of the current experimental constraints 
are satisfied. We consider constraints from DAMA/Nal- 
96 @, DAMA/NaI-03 0, CDMS-I @, CDMS-II 0, 
CRESST-I 0, and CRESST-II 0. The exposures, ef- 
ficiencies, thresholds, and constraints we use are listed in 
Table 1. 

We compute the expected number of recoil events with 
recoil energy in the range (Ei,E 2 ) as 

i 1 

The sum runs over the nuclear species in the detector, 
and Si = MiT.iei(E) is the effective exposure of species i 
(Tj being the time the mass Mi is exposed to the signal, 
and €i(E) being the efficiency). Moreover, dRi/dE is the 



expected recoil rate per unit detector mass and unit time, 

dRj = p<Ji\Fj(E)\ 2 j f{y,t) ^ 

Here Mj is the nuclear mass, fa — mMj/(m + Mj) is the 
reduced WIMP-nucleus mass, p is the local halo WIMP 
density, Fi(E) is a nuclear form factor, Oi is the WIMP- 
nucleus cross section, and /(v, t) is the WIMP velocity 
distribution in the reference frame of the detector. 

In this analysis, we assume that the WIMP-nucleus in- 
teraction is spin-independent and scales with the square 
of the nucleus atomic number Aj as Oi = apAf(pi/p p ) 2 . 
Here a p is the WIMP-proton cross section. For the nu- 
clear form factor we use the conventional Hclmi form, 
Fi(E) = 3e~ q s / 2 [s\n{qr) — qr cos(qr)]/(qr) 3 , with s = 
1 fm, r = VR 2 - 5s 2 , R = 1.2A 1 / 3 fm, q = v / 2M l E. 

To find a compatible region, we vary the WIMP mass, 
the stream density, and the stream velocity (allowing also 
for the absence of the stream). For each choice of these 
variables we find the WIMP-proton cross section a p that 
produces the desired modulation amplitude in the 2-4 
keVee bin of DAMA/NaI-03 within two sigma (the re- 
sulting <7 P are shown in Fig. 2; we tried other ways to 
fix cr p using also the 2-5, and 2-6 keVee bins and ob- 
tained similar results). With this WIMP-proton cross 
section, we evaluate the expected number of events in all 
of the other experiments. We thus determine if the pa- 
rameters we choose are compatible with the experimental 
constraints we impose. 

Fig. 1 shows the compatible region in heliocentric 
stream speed versus WIMP mass resulting from our 
study cut arbitrarily to u str < 1200 km/s. The region is 
bounded from above by the constraints from CDMS-I and 
II, and from below by the constraint from DAMA/NaI-96 
and CRESST-I, up to masses of about 6 GeV. The ver- 
tical band around m = 6 to 9 GeV shows that there the 
stream speed can be arbitrarily small. Indeed, in that 
mass range, the DAMA modulation can be reproduced 
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FIG. 2: Region of WIMP-proton cross section a p resulting 
from our study vs. WIMP mass m (cross-hatched region: 
without stream; hatched region: with stream and velocity 
cuts as shown). Also limits from DAMA/NaI-96, CRESST-I 
and II, and CDMS-I and II (see Table 1). 

with just the conventional Maxwellian distribution, with- 
out an additional stream. 

The hatched region in Fig. 2 shows the values of the 
WIMP-proton cross section o~ p resulting from our study 
that are compatible with the DAMA/NaI-03 modula- 
tion, as a function of the WIMP mass. The region was 
cut arbitrarily at i> str < 1200 km/s, as in Fig 1. Also 
shown is where the cut would be for a Galactic stream, 
i> s tr < 840 km/s. The cross-hatched region is the case of a 
conventional Maxwellian halo without a stream. The up- 
per bounds on a p are given by CRESST-I, CDMS-I and 
II, and DAMA/NaI-96. The values we find are included 
in the region presented as allowed by DAMA in Fig. 28 



of Ref. . Li ght neutralinos as WIMPs with masses as 
low as 2 GeV |15) or, with updated bounds, 6 GeV 
have been considered, but their cross sections are about 
one order of magnitude smaller than those needed here. 

In conclusion, we have pointed out that for light dark 
matter particles a signal could be observed by DAMA 
through its Na component. This signal would be below 
threshold for Ge in CDMS and EDELWEISS. This pos- 
sibility can be tested with a few months of Si data in 
CDMS-II, and future O data in CRESST. 

For WIMPs with spin independent interactions, we 
have presented two examples of dark matter velocity dis- 
tributions that give the annual modulation observed by 
DAMA but satisfy all other constraints from dark matter 
searches. The first is a conventional Maxwellian distri- 
bution with a WIMP mass around 6 to 9 GeV. Surpris- 
ingly, this simple possibility remains open. Our second 
example is the former distribution superposed to a dark 
matter stream coming from the general direction of the 
Galactic rotation (perhaps associated with extragalactic 
dark matter, but not the Sagittarius stream). 

For the sake of illustration, we have assumed a partic- 
ular density of the stream (0.03 of the local halo density) 
and a particular incoming direction (halfway between the 
plane of the Earth's orbit and the direction of the Sun's 
velocity in the Galaxy). Our main result are the allowed 
regions presented in Figs. 1 and 2. 

Clearly, the effect of the stream would be larger 
(smaller) for larger (smaller) stream densities and for in- 
coming directions closer to (further from) the plane of 
the Earth's orbit. For simplicity, we have illustrated our 
idea only for the case of WIMPs with spin-independent 
interactions. Other kinds of particles and interactions, 
or halo velocities distributions more complicated than a 
conventional Maxwellian distribution, may extend the al- 
lowed regions of parameters. 
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